General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



ACUITY FORM 602 


X-641-70-317 

PREPRINT 


nasa m x- 653 o f 

KINETIC THEORY OF PLASMA IN 
A MAGNETIC FIELD 


JOHN C. PRICE 


AUGUST 1970 


GODDARD SPACE FLIGHT CENTER 


GREENBELT, MARYLAND 

* N 7 () - n fi 5 flfi 


— — - — ^ tj \J . — 

(ACCESSION NUMBER) 

Jy 

— 1 JJ U 

(TH^I) 

(NASAtR OR TMX OR AD NUMBER) 

(CODE) 

2 % 
(CATEGORY! 
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John C. Price 
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Goddard Space Flight Center 
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Green be It, Maryland 20771 

ABSTRACT 

The Krylov-Bogoluibov transformation is used to simplify the 
guiding center equations for a charged particle in electric and 
magnetic fields. The resulting transformation of variables is then 
applied to the Vlasov distribution function, yielding a magnetic 
Vlasov equation which describes the low frequency behavior of a 
system with no statistical effects. The equation is generalized 
to include effects of high and low frequency fluctuations by a 
procedure developed by Klimontovich and Dupree. A consistent 
treatment of the conservation laws and Maxwell's equations is given 
to complete the kinetic description. 

I. INTRODUCTION 

Many experiments dealing with plasmas in a magnetic field take place in 
the difficult regime where the collision mean free path of a particle is 
comparable to, or longer than, the characteristic scale lengths (gradient, 
curvature, etc.) of the system. For these situations it is desirable to have 
a theory which contains both a partial solution of the equations of motion 
of a single particle and the corrections due to statistical fluctuations, in- 
cluding collisions due to particle discreteness. 

The standard method of solving the particle equations for motion for 



r(-i) j'V'tt) introduces complications when a strong magnetic field is 
present, for the rapid oscillation about the field line produce a lengthy 
expression for the orbit. This in turn causes great complexity in the stat- 
istical theory of many particles, even for the most simple geometry . 1 

Recently, Wilson 2 has carried the analysis of guiding center motion 
through order £ , and has written the conservation equation for guiding 

centers. This work is conceptually similar to his, and to that of Hastie 
and Taylor 3 , but with considerable simplification in detail and with a 
development which connects with earlier work in statistical kinetic theory . 4 
In addition, Wilson 5 has reviewed much of the guiding center work to date. 

In Section II by using the Krylov-Boglouibov expansion technique and 
several modifications of the definition of initial guiding center variables 
we produce economical expressions for single particle motion. The Krylov- 
Bogoluibov method has the virtue that the new variables are defined to be 
equal on the average (over the phase angle about the magnetic field) to the 
old variables, so that an intuitive identification is possible. 

In Section III this change of variables is applied to the Vlasov 
equation, yielding a magnetic Vlasov equation which is valid for low frequency 
disturbances if statistical effects may be ignored. Because the definition 
of the distribution function is precise we are able to identify the charge 
and current for Maxwell’s equations without the need to carry auxiliary 
moment equations as in Wilson's work.® 

Section IV develops the effect of statistical fluctuations, including 
the definition of the statistical transformation of variables, and the 
equation of motion for fluctuations. Although the method described 
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earlier 4 is now applicable, we defer evaluation of the kinetic equation 
until a specific problem makes further approximation possible. 

II. SINGLE PARTICLE MOTION 

The motion of a particle of charge ^ and mass r*\ is described by the 
equations 

r - v* 

& L 


where the fields' are in Gausian units and £ L . in the usual way 7 

we wish to expand in € , treating the revolution about the magnetic field 
lines as the lcrwes-t order effect. Throughout the paper we work' Only to first 
order in £ , so that we may omit the explicit rescaling of time which limits 
the guiding center theory to frequencies far below the cyclotron 1 frequency. 8 
We define the: guiding center variables ( ft , 0 ) 


r - R + ^[Nfe/Occ^ -m 


ecw) 


'V h = v L (f?A) 

■v^- (v'f- -vf 


V” ‘V- l( 



are 


parallel 


Here L , M , A/ are orthogonal unit vectors with i. ( {-) 

to the magnetic field at R , and and A / chosen in any convenient way 
with A/ - L* M . The vector U satisfying (J 0 remains to be 

determined. 

These variables satisfy the equations 


• ("/}-♦- L 
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where for a quantity £ (R,-f ) j + R V C and 

A -- c e (*,*> - u ' ems) - *1%*^%] 
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1 The parallel component of B is assumed to be of order £ , i.e., 

^ We now choose U such that the perpendicular component of A equals Q . 

f ; 

§> Thus - U 0 ■* € (A, - ■ • . where 

r 

c £(£,+') 'L(RJ) 

dfe,+) 




€ 


LKJ) 
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f Also we Taylor expand about using the relation K-R- -~(bl c t'*Q -Afs /'►,£) 

| in order to eliminate ^ . For the accuracy we require it is sufficient- 

f to keep angle independent terms through order € , and periodic terms through 

BT 
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Vj. - -j[-ct--pxr* (u+V'O iB ] + ^ [(*f«'M) viu+v,L) ccsJ4 + 
Ain) : q(u*v„ 

Q - - | ~~ r ^ tjffi + ~z L v* /tt+ y ’» L )] ~ (bIoo%Q -M $>*&)’ vB 

j|(am^ / i/aj) : v/iA+m, *■) Co\>2Q - (M/i-AJA))-’ q/Ui v>, l) s>* 2&J 

+ 6(s) 


(3) 


The order € part of £ is lengthy but will not be needed. Also, the 
angle independent terms of Vi. are equal to (yi /ag) ^ 


We now use the Krylov-Bogoluibov method 1 °> 11 to define new variables 

0,<, 0 ) -* ( P, *H 2, $) 

R- P 
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where all quantities on the right are located at P , e.g . ^ (P,+). These 
variables satisfy the I independent equations 


P: U*WL 
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(c) 
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These transformations of variables lead to considerable simplification in 
the kinetic description of a plasma. 
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III. MAGNETIC VLASOV EQUATION 
If the distribution function 
volume "y satisfies the equation 


normalized to 





( 6 ) 


where the fields satisfy the restriction of Section II, then the trans- 
formations (l), (4) produce a different functional form. Because of 
equation 5 C and the fact that ^ is positive we make the further substitution 
°yy\= z* Aws) define 

( 7 ) 

By chain rule differentation F satisfies the equation 

PaF J '0 

ft 

(S) 
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where the coefficient? rf the derivatives are given by 5, valid to order € . 

Typically we assume p varies slowly in tine, so that the $ term may be 
dropped. 

Equation 8 plus Maxwell's equations describe the behavior of a system 
where the initial conditions are known and for times short enough so that 
statistical fluctuations are not important. In order to complete Maxwell's 
equations wc require the charge density and current, which may be developed 
from 7 in the following way. 

We carry out the velocity integration 12 in the original guiding center 
variables Q 




^Current density 




r Hrfh 
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Her*“ "F, is the system average density. 


r 2 
Now we expand p about the variables f*. ^ Q 


Thus 
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where the order £ quantities (P'r), (%l- V.,) , etc., are 

given by the transformation equations 1 aad 4 and the definitions of and 
hf\ , and F ( (,V ll/ m / 0 y ^ ) is simply a relabeling of f ( 

$ t) . Since p is periodic in Q- it may be 

expanded in Fourier series; the requirement that fields vary slowly in time 
means that the contribution to K and T from f- H t must vanish. 
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in order 


Finally we rename the dummy variables of integration ^7 ,^4 , S& 
to produce a notation consistent with 8. This is not a further change of 
variables, but simply a relabeling. We have 


K(r, *) = 


<P- 
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7 6o +) = fir 2 * | <r3/f[U+1/L * (HiN-Al vM) 


-f 


6 


L*v 6 - 





(11) 


where in both cases r is the angle independent (average) part of the total 
distribution. Maxwell's equations in r . 4r are given by 

r , / 30 

y* * ~ c. 


( 12 ) 
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IV. STATISTICAL THEORY 
A. Conservation Lavs 

In reference 4 a procedure was developed for treating the 
Klimontovich-Dupree hierarchy of equations. In order to correct a defect 
in the conservation properties of that work we write the equation for the 
single particle distribution of species , neglecting electromagnetic 
effects in the fluctuating fields, i.e., $8 Z Q . 




^ o r* <6 )1 
c J Jy 



( 13 ) 


Here the brackets ) mean a statistical or ensemble average. Although 
we shall change variables in order to eliminate rapid phase dependence, 
equation 13 as written is convenient for developing conservation of energy and 
momentum in the system. In this paper the state of the system is described by 
the one particle distribution, and by the distribution of electrostatic energy 
in local modes, i.e., fluctuations in which the perturbed electric field may be 

approximated locally by ^ /<F£fc) €*p (*■ K‘ * ~ K 

with - O-ic \ . Thus we assume knowledge of -Ac c t) 

and ^ (r ^ ~ ^ t snd. for a11 

time on any physical boundaries. 
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Multiplying the right side of 13 by V and ^ ^ 7 

integrating by parts, and using Maxwell f s equations in the longitudinal 
approximation yields 14 




^ V <f 7T / 






( 14 ) 


We adopt these equations for the determination of D-p. and T*. 
instead of determining them from the dielectric function. In the homogeneous 
field free case these equations reduce to the imaginary and real parts of the 
usual dielectric function . 15 These relations guarantee momentum and energy 
conser/ation in the system, where the energy is given by J^/vr -f 4 


£ 

* 




B. Statistical Variables 

In this section we utilize the transformations of Section II to develop 
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the statistical theory of plasma in a magnetic field. Because of fluctuations 
in the electric field we define average variables before following the pro- 
cedure of Section III. 

We assume that instruments which measure electric fields are able to 

i 

disregard or average over high frequency (00 > Cj^ ) or short wavelength 
fluctuations, but that they respond to fluctuations which meet the require- 
ments of the guiding center theory. Then we may define ensemble average 
variables by splitting the electric field (measured) = £ 

where is the fluctuating guiding center field. We define 

R- R - X e f e 

— r-r 3 " 

'Vu ' ^ ® 6 te 

•Vi. - T*. - 

e-- e- 


( 15 ) 


We shall disregard terms which lead to results of order <S C S ^ 
in the final equations. 

As in Section II we calculate (r)^ (■$■„) , etc., and 

use the Krylov- Bogoluibov transformation to define new variables { p , tj , jr, ^ ) 
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These satisfy the equations of motion 


p - P(f'%*A) 
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( 17 ) 
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where the notation indicates simply the relabeling of the right side of 
equation 5. However, the electric field is the sum of the average field and 
fluctuations which satisfy the guiding center restrictions. We ensemble 
average 17 to find 


</>y -- U-* i L * 

<*> ' ^ ■ «■£ - * U'L (l * * f >/l Vi «t)\ 
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where L\ -- + $ U, 
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By subtracting (l8) from (17) and defining , etc.. 
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ft- 

e 



{ 7 C vB 

1B L 


Jf -- 7 



(19) 


C . Statistical Equations 

We now use the method developed by Klimentovich and Dupree 17 and the 
transformations developed in this paper to write statistical equations de- 
scribing the behavior of a plasma in a magnetic field. The exact one particle 
distribution 


M 


f( r s V,*)- ~ E 'b(r-Yjt)') J/Vvft m) 
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satisfies the equation 


If 




Here £ is the exact electric field, including, for example, particle 
discreteness effects and high frequency collective fields 
Defining £" z *■ J 4 ^ hi and using the change of 

YiBU.i) wetoTO 


variables 15, 1 6, followed by ^ - * 

~ 




[^T enzr — h 
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We define the ensemble average to have no ^ dependence 


F (w-i'i'y = -f(p,i y/ a y {) 


and take the average of 20, assuming high frequency and low frequency modes 
are separable (e.g., by Fourier analysis) so that F~ f - 

l '<r> ’ * wi,] f ■ - < si. »#„> - 






/cfr> 4 -j>) <? <£& 


6 & 9^ 
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where 





In the notation of reference 4 we subtract to find the equations for 
the fluctuations. 




li . <i> , * <» ? 7 * <*> %])*«> -- - ^ 

* <?*J L 1 ^ e 5* 





Maxwell’s equation uay be developed as in Section III. For the 


longitudinal fluctuations we require only the charge density, 


V IK> - 


We may now use methods developed earlier 4 ' 16 to solve equations 22 end 23 
by integration along the characteristics, and insert the results into 21, The 
low frequency terms represent a generalized form of Dupree’s vor£? while the 
high frequency terms generalize the result of Rostoker 1 . Because of the 
great length of the resulting equation we omit the expression until 
application to a specific problem makes further approximation possible. 
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